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Weak Decay Form Factors of Heavy Pseudoscalar 
Mesons within a Light-Front Constituent Quark 



The transition form factors describing the semileptonic decays of heavy pseudoscalar 
mesons are investigated within a relativistic constituent quark model formulated on the 
light-front. For the first time, the form factors are calculated in the whole accessible 
kinematical region, adopting meson wave functions derived from various effective qq 
interactions able to describe the meson mass spectra. It is shown that the decay rates 
of the B — > Diui, D — > Kivg and D — > nii'i weak decays are mainly governed by the 
effects of the confinement scale, whereas the B nii'i decay rate is sensitive to the 
high-momentum components generated in the meson wave functions by the effective 
one-gluon-exchange interaction. Our results are consistent with available experimental 
data and predictions of lattice QCD calculations and QCD sum rules. 
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The investigation of semileptonic heavy-quark decays is an important test of our 
understanding of weak and strong interactions. As is well known, the semileptonic decay 
amplitude factorizes into the product of the leptonic and hadronic V—A currents. The matrix 
elements of the latter contain relevant pieces of information about the strong forces which 
bind quarks and gluons into hadrons, whereas the leptonic part depends on the Cabibbo- 
Kobayashi-Maskawa {CKM) mixing parameters, Vq-^q^, which are fundamental quantities of 
the Standard Model. In order to calculate the decay rates and, hence, to extract the CKM 
parameters from the experiments, it is necessary to know the hadron form factors in the 
whole region of accessible values of the squared four- momentum transfer, < < q^ax- 
should be pointed out that till now the predictions obtained within various non-perturbative 
approaches, like, e.g., the QCD sum rule technique and the constituent quark (CQ) model, 
do not cover the full range of g^. In particular, as far as CQ models are concerned, the weak 
hadron form factors are usually calculated at a fixed value of (appropriate for the specific 
CQ model) and then extrapolated to the whole range of q^. As a matter of fact, in the 
original WSB approach the form factors are calculated at g^ = and then extrapolated 
to g2 > using a monopole ansatz 1/(1 — g^/Mp^^g), with Mpoie being the mass of the lowest- 
lying vector meson in the given channel. Within the ISGW model the form factors are 
calculated at the point of zero recoil (i.e., g^ = q'^ax) then extrapolated to g^ = using 
an exponential-like ansatz. The CQ model of Ref. P] uses the light-front (LF) formalism 
to compute the form factors for space-like values of g^; then, the extrapolation to the time- 
like region is performed using a two-parameter formula, which reproduces the values of the 
form factors and their first two derivatives at g^ = 0. Finally, in Ref. |^ the form factors 
have been calculated in the space-like region using a dispersion integral approach and, then, 
analytically continuated in the time-like one. 

In Ref. the g^ dependence of various heavy-to-heavy and heavy-to-light form 
factors among pseudoscalar mesons has been calculated in the whole accessible kinematical 
region. The calculations have been based on a LF CQ model, adopting a gaussian-like ansatz 
for the meson wave functions, and carried out in a reference frame where the momentum 
transfer is purely longitudinal, i.e. in a frame appropriate for time-like values of g^. It has 
been shown that the time-like LF result, obtained using the matrix element of the "good" 
component of the weak vector current = J^ + J^, coincides with the contribution of the so- 
called spectator pole of the Feynman triangle diagram. The remaining part of the Feynman 
diagram, the so-called Z-graph, corresponds to the contribution of quark-pair creation from 
the vacuum. The Z-graph contribution is expected to be marginal for heavy-to-heavy decays 
(see PI) and to become more important near the zero- recoil point for heavy-to-light decays 
(see 0). In this letter we use the approach of Ref. P, but, as far as the wave functions 
of the initial and final pseudoscalar mesons are concerned, we adopt the eigenfunctions of 
LF mass operators, constructed using various effective gg interactions able to reproduce the 
meson mass spectra. In this way a link between heavy-meson weak decay properties and the 
"spectroscopic" constituent quark model is explicitly constructed. It is shown that the decay 
rates for the B —>■ Diu^, D Kiui and D niui decays are mainly governed by the effects 
of the confinement scale, whereas the B —>■ niu^ decay rate and, hence, the extraction of \Vub\ 
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are sensitive to the high-momentum components generated in the meson wave functions by 
the effective one-gluon-exchange (OGE) interaction. Our results obtained both for the form 
factors and the decay rates are compared with available experimental data and predictions 
of various non-perturbative approaches, like, e.g., lattice QCD simulations and QCD sum 
rules. 

Here and below we denote by Pi (P2) and Mi (M2) the four-momenta and masses of 
the parent (daughter) mesons, respectively. The momenta and masses of the corresponding 
active quarks will be denoted as pi (^2) and mi (7/12), whereas k (m) is the momentum 
(mass) of the light spectator (anti) quark. The four-momentum transfer is g = Pi — P2 and 
the fraction y is defined a.s y = P2 / Pi ^ where P^ = P^ + P^ with the 3-axis defining the 
spin quantization axis. In a frame where the three-momentum transfer is purely longitudinal 
(i.e, = 0), it can be easily verified that: = (1 — y){Mf — M^/y). Thus, one gets: 
1/1,2 = Cijl i V'?^ ~ 1); where C = M2/M1 and rj = Ui ■ U2, with Ui and U2 being the four- 
velocities of the parent and daughter mesons, respectively. The two signs for ?/i 2 correspond 
to whether the three-momentum transfer is in the positive or negative direction of the 3-axis. 
At the point of zero recoil, = g^^^ = (Mi — M2Y, one has yi = y2 = while yi = I and 
1/2 = at g2 = 0. 

The semileptonic decay of a pseudoscalar meson Qiq into another pseudoscalar meson 
Q2q is governed by the weak vector current only. The transition form factors h±{q^) and 
f±{q^) are defined as 



< P2\Q2lMPi > = ^MiM2 [h+{q') {Ui + U2), + h_{q') {Ui - U2)^ 

= U{Q'){Pi + P2), + f-{Q'){Pi-P2), (1) 
Neglecting lepton masses, the inclusive semileptonic decay rate F is given by 

r = |i'«.*P^C'f~ A, (r-ip /?(.,) (2) 

where r/max = (C + 1/0/2 and Vq^q^ is the relevant CKM matrix element. The form factors 
/±(Q'^) and h±{q^) can be evaluated using only the "good" component of the current J"*" = 
Q^l^Qi- The matrix elements of depend upon both g^ and y, whereas the form factors 
/± (^±) are independent of y. In order to invert Eq. (p, two matrix elements J~^{q'^, y) =< 
P2\Q2'l~^Qi\Pi >, corresponding to y = yi and y = y2, should be calculated, obtaining: 
J+(g2,i/,) = 2P+ [/+(g2) (l + ^,) + /4g2) ^ith I = 1,2. Putting H{q^,y,) = 

J^{(f,yt)/'^Pi one has 

^^^^2^ _ i.y2 T l)H{q^, yi) - {yi T l)H{q^, ^2) 



Following Ref. 0], the LF ground state of the parent meson Qiq can be written as 
dk± dx 



Pi >= E / , , A(^ (X, kr, XX) ^ QU^O > (4) 
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where Q\ = Qi{pi, A) and = q^{k, A) are the creation operators of the heavy active quark 
and the hght spectator antiquark, with hehcities A and A, respectively. The LF momenta 
are taken as: k = [k^, k±) and pi = Pi — k = {Pi' — k~^, —k±), where k± is the relative 
transverse momentum. The spectator and the active quarks carry the fractions x and 1 — x of 
the plus component of the momentum of the meson, respectively. In Eq. (^) i?oQ^(x, k±] AA) 
is the momentum-dependent spin part of the meson wave function, arising from the Melosh 
rotations of the CQ spins (see 0), and Ai{x, k^) = Mio[l — (mf — m^)^/Mj^Q] / 4x{l—x), with 
M^Q = {ml + k'^)/{l — x) + {rn? + k\)/x being the squared free mass. Finally, k"^ = k'^ + fcg, 
where k^ = {x — l/2)Mio + {ml — m^)/2Mio. The state vector of the daughter meson has 
a form analogous to with the obvious replacement 1^2. The states \Pi > {i = 1,2) 
are normalized as: < Pl\Pi >= (27r)^ 2Pj^ ^{Pt ~ P't)^ so that the normalization condition 
for Wi{k'^) is /q°° dk k'^wf{k^) = 1. Putting x' = x/y for the LF momentum fraction of the 
spectator quark in the final meson, the matrix elements H{q^, y) can be cast in the following 
form 

H{q\y) = I'dx I dk^ ^Ai{x,kl) Mx',kl) 



[m{l — x) + xmi] [m{l — x') + x'm2] + kj 



^[m{l — x) + xmi]2 + k'j_ ^[m{l — x') + a;'m2]^ + k\ 



(5) 



The radial wave function Wi{k'^) can be chosen to be the eigenfunction of an equal- 
time effective Qq Hamiltonian (see for more details 0). For the latter we will consider 
two choices: the first one corresponds to the relativized Godfrey-Isgur {GI) Hamiltonian 
|§] and the second one is the non-relativistic {NR) model of Ref. ||^. We want to point 
out that both potential models nicely reproduce the mass spectra of light as well as heavy 
mesons. In particular, the interaction terms contain an effective OGE interaction, composed 
by a spin- dependent part, responsible for the hyperfine mass splitting in light mesons, and 
a spin-independent term, responsible for the hydrogen-like pattern of heavy-meson (botto- 
nium) states. For comparison, we consider also a Gaussian-like ansatz, with the harmonic 
oscillator {HO) length taken from the updated version of the ISGW model |]10|- In what 
follows, the three choices will be referred to as the GI, NR and HO cases, respectively. The 
CQ momentum distribution corresponding to the wave functions w^'^^\ w^'^^^ and w^^^\ 
obtained for the D-meson, is shown in Fig. 1. It can be seen that: i) w^^'~'\k'^) is quite 
similar to the wave function obtained by retaining only the confining part of the GI inter- 
action, i.e. it takes into account only the effects of the confinement scale; ii) both w^^^^ and 
yji!^R) exhibit high-momentum components generated by the effective OGE term; iii) the 
high-momentum tail is higher for the GI interaction. Similar results (cf. 0) hold as well in 
case of the other mesons (vr, K and B) considered in this letter. 

The results obtained for the form factor f+{q^) (Eqs. (^ and (|^)) and the decay rate 
F (Eq. (0)) of the B —>■ Dlu^, B —>■ iriui, D Kiui and D —>■ Triui decays, are reported in 
Tables 1-3 and Figs. 2-3, and compared with results from lattice QCD calculations, QGD 



sum rules and quark models. In our calculations the PDG |Tl| values of the meson masses 



have been adopted for all the three models of the meson wave function. 
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Decay B Diui. From Table 1 it can be seen that our results for f^~'^{q^ = 0), 
f+^^ilmax) the decay rate T{B — >• Diui) depend only slightly upon the choice of 
the meson wave function, i.e. these quantities are mainly governed by the effects of the 
confinement scale. Our calculated rates are in agreement with existing CQ model predictions, 
which typically range from ~ 7 to ~ 13 ps~^. By combining the branching ratio Br{B^ 
D-i+Ui) = (1.9 ± 0.5)% with the world average of the B° lifetime, tbo = 1.50 ± 0.11 ps 
11|, the measured rate is r(i?° —>■ D~l^vi) = (1.27 ± 0.35) ■ 10^° s~^. From our predicted 



rates (see Table 1), one gets: \Vbc\ = 0.034 ± 0.005 (GI) and |Hc| = 0.036 ± 0.005 {NR 
and HO). These predictions, which are expected to be modified only slightly by radiative 
corrections (see |12|)! are in agreement with the updated "experimental" determinations of 
\Vbc\ obtained from exclusive and inclusive semileptonic decays of the B-meson, viz. \Vbc\exci = 
0.0373 ± 0.0045e^p ± 0.0065*^ and |Hc|mcZ = 0.0398 ± O.OOOSe^p ± O.OOAOth 0. At the point 
of zero recoil the form factor h+{ri = 1) results to be: 0.960 (GI), 0.971 (NR), 0.965 (HO), 
in nice agreement with the updated QCD sum rule estimate, 1 + Si/m^ = 0.945 ± 0.025 
T2| . Moreover, the values of the slope, p^, and the curvature, c, of the form factor h^{ri) 



at the point of zero recoil, i.e. h^{ri) ^ /i+(l){l — fP{ri — 1) + c(?7 — 1)^) + 0[{ri — 1)^]}}, 
are: = 0.96 (G/),1.22 (iVi?), 1.30 {HO) and c = 0.71 (G/),0.96 (A^i?),1.20 (HO). The 
results corresponding to the GI wave function compare very favourably with the recent 



model-independent improved bounds obtained in [Q. The full dependence of the form 
factor f^~'^{q'^) is shown in Fig. 2(a). It can be seen that: i) in the whole range of 
the form factor f^^^{q^) is affected only slightly by the high- momentum tail of the B- and 
Z?-meson wave functions; ii) the monopole approximation, with Mpoie = 6.25 GeV [|T|, yields 
a slightly flatter q^ dependence. 

Decay B nii>£. The dependence of the form factor f^~'^{q^) is shown in Fig. 
2(b). Near the zero-recoil point our result does not exhibit any pole dominance and deviates 



from recent results of lattice QCD simulations |T5|. Such a discrepancy can be ascribed 
to the effects related to the admixture of the \B*'k > component in the i?-meson wave 
function 0, which, in the LF language, corresponds to the contribution of the Z-graph 
(quark-pair creation diagram). We want to point out that the decay rate T{B — >• niue) is 
not expected to be sharply affected by the Z-graph, because in Eq. (0) the contribution 
arising from the region near the zero-recoil point is kinematically suppressed. From Table 
1 it can clearly be seen that the form factor f^^'^{q^ = 0) and, hence, the decay rate 
r{B — s> niui) are sensitive to the high- momentum components generated in the meson 
wave functions by the effective OGE interaction, particularly in case of the GI model. 
This feature is related to the huge value of the maximum recoil achieved in the B — > Tciui 
transition. Moreover, our LF results are in overall agreement with the ones obtained from 



recent lattice QCD simulations [|l^, |T^, and only slightly lower than the results of the 



analysis performed in Ref. [0] using the LF QCD sum rule (see Tables 2-3). Recently, 
the CLEO collaboration |20| has reported the first signal for exclusive semileptonic decays 
of the B meson into charmless final states, in particular for the decay channel B Tciui. 
Adopting the models of Refs. |]l| and [0] to reconstruct the efficencies, the "model-dependent" 
branching ratios are: Br{B iriue) = (1.63±0.57)-10~^ and (1.34±0.45)-10~^, respectively. 
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Combining the average of these results with the world average value of the 5° lifetime, one 
gets: r{B^ Ti~C.^Vi/) = (0.99 ± 0.25) • 10~^ ps~^. From our predicted rates (see Table 1) 
one has: \Vbu\ = 0.0025 ±0.0003 (G/), 0.0030 ±0.0004 (NR) and 0.0032 ±0.0004 (HO). This 
means that the extraction of \Vub\ is sensitive to the effects due to the high- momentum tail 
of the meson wave functions. We stress that this conclusion is not expected to be sharply 
modified by the effects related to the quark-pair creation from the vacuum. Using the results 
previously obtained for \Vbc\, one has: \Vbu/Va\ = 0.074 ± 0.014 (G/), 0.083 ± 0.016 (NR) 
and 0.088 ± 0.016 (HO). These results are in accord with the value of iHu/Kbl derived 
from measurements of the end-point region of the lepton spectrum in inclusive semileptonic 
decays |2l], H, viz. = 0.08 ± O.Ole^^p ± 0.02^;,. 



Decay D Kiui. The behaviour of the form factor f^^^{q'^) is shown in Fig. 
3(a). It can be seen that our form factor is mainly governed by the effects of the confinement 
scale and it is consistent both with the assumption of pole dominance and with the results 
of the lattice QCD calculations of Refs. [1^, Moreover, our results for f^^^{0) and 
f+^^ilmax) (sse Table 1) compare favourably with the recent experimental results p4| 



f^-'^iO) = 0.77 ± 0.01 ± 0.04 and f+^^{qiax) = 1-42 ± 0.25, obtained from the measured 
lepton spectrum assuming pole dominance. By combining the experimental data on the 
branching ratio, Br{D^ — > K~e~^i'f.) = 3.68 ± 0.21% |TT[, with the accurate value of the 
lifetime, t{D°) = 0.415 ± 0.004 ps one has: T{D° i^"e+z/e) = 0.089 ± 0.005 ps-\ 
From our predicted rates (see Table 1) one obtains: \Vcs\ = 0.88 ± 0.03 (GI), 0.90 ± 0.03 
{NR) and 0.89 ± 0.03 (HO). The constraint of unitarity of the CKM matrix with three 
generations of leptons yields a (~ 10%) higher value, viz. \Vcs\ = 0.974 [O . 



Decay D — ^ niui. This is the only heavy-to-light decay where an extensive com- 
parison with experiment is possible. Our results for the form factor f^^^{q^) are shown 
in Fig. 3(b). It can be seen that its behaviour at low agrees with the prediction of 
the LF QCD sum rule of Ref. [l^, whereas near q^ = g^^^ it deviates from the pole ap- 



proximation. The same discussion, already done about the relevance of the Z-graph in the 
B —>■ IT decay, applies as well to the D —>■ n transition. From Tables 2-3 it can be seen 
that our values for f^^'"{q^ = 0) and r{D niiye) are in nice agreement both with the 
experimental data [|ll|] and the results obtained within various non-perturbative approaches. 
Assuming \Vcd\ = 0.221 ± 0.003 |TT[] from the unitarity of the CKM matrix, our predic- 
tions for the semileptonic decay rate are: r{D^ TT~e^Ve) = 8.2 ± 0.2 (CI), 7.6 ± 0.2 
{NR) and 7.8 ± 0.2 (HO), in units 10~^ ps~^. These results compares favourably with the 
experimental value (9.412^9) ■ 10~^ ps^^ [11] and the recent LF QCD sum rule prediction 
(7.6±0.2)-10~^ ps~^ |l^. The ratio of the branching fraction of the Cabibbo suppressed decay 
D Tilvf to that of the Cabibbo favoured decay D Kiui has been recently determined by 



the CLEO collaboration by measuring both charged p5| and neutral p6[ D-meson decays. 
Assuming the pole dominance for the form factors, the following values have been obtained: 
I/+ ^^(0)7/^*^^(0) I = 1.29 ± 0.21 ± 0.11 H and 1.01 ± 0.20 ± 0.07 [||]. With respect to 
the average of these values we predict a slightly lower ratio, viz. |/+~"^(0)//:P^-^(0)| = 0.91 
(GI), 0.88 (NR) and 0.87 (HO). The value obtained in Ref. ^ is 0.71 and other theoretical 
predictions typically range from 0.7 to 1.4. 
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In conclusion, adopting a light-front constituent quark model we have investigated 
the transition form factors, which govern the heavy-to-heavy and heavy-to-light semileptonic 
weak decays between pseudoscalar mesons. For the first time, the form factors have been 
calculated in the whole kinematical region accessible in semileptonic decays, adopting meson 
wave functions derived from various effective qq interactions able to describe the meson 
mass spectra. It has been shown that the decay rates for the B Diu^, D Kiu^ 
and D —>■ ttIui weak decays are mainly governed by the effects of the confinement scale, 
whereas the B — > iriui decay rate and, hence, the extraction of \Vub\ are sensitive to the 
high-momentum components generated in the meson wave functions by the effective one- 
gluon-exchange interaction. Our results both for the form factor and the decay rates are 
consistent with available experimental data and predictions of lattice QCD calculations and 
QCD sum rules. We want to stress that an estimate of the contribution of quark-pair 
creation from the vacuum is mandatory, particularly in case of the B niui and D — » tt£i>i 
transitions. 

We are very grateful to K.A. Ter-Martirosyan for valuable discussions. Two of us 
(I.L.G. and I.M.N.) acknowledge the financial support of the INTAS grant No 93-0079. This 
work was done in part under the RFFR grant, Ref. No. 95-02-04808a. 
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Table Captions 



Table 1. The form factor f+{q^) (Eqs. (j^) and (^) for various semileptonic decays, eval- 
uated at = and = q"^^^ using the GI, NR and HO wave functions (see text), 
and the corresponding decay rate T (Eq. @), calculated in units ps~^ assuming 

I^QiQ2l = 1- 

Table 2. The form factor f+{q^) of the B -nivn and D niui transitions evaluated 
at = within different approaches. The label Exp means the experimental result 
(assuming pole dominance), while the labels SR and LAT correspond to QCD sum 
rule and lattice QCD calculations, respectively. 

Table 3. The same as in Table 2, but for the decay rate F of the B niui and D niug 
transitions, calculated in units lO^^s"^ and 10^^s~\ respectively, assuming IVq^QjI = 1. 
The experimental result quoted for the D — > niui transition has been taken from [jll| 
assuming \ Vcd\ = 0.221. The results from the lattice QCD simulations of Refs. [T^, |16| 
are obtained assuming pole dominance. 



Figure Captions 

Fig. 1. The CQ momentum distribution \w{k'^)\'^ versus the internal momentum k in the 
D-meson. The solid, dashed and dotted lines correspond to the CI, NR and HO 
cases, respectively (see text). The dot-dashed line is the result obtained using the 
eigenfunction of the Hamiltonian of Ref. [§] in which only the confining part of the 
interaction term is retained. 

Fig. 2. The form factor f+{q^) of the B — >■ D^v^ (a) and B —>■ Tiiui (b) transitions. The 
solid, dashed and dotted lines correspond to the results of our LF calculations (Eq. 
(H)), obtained using the CI, NR and HO wave functions, respectively. The dot- 
dashed lines are the monopole approximation f+{q^) = /+(0)/(l — q'^/M^^iJ with 
Mpoie = 6.25 CeV (a) and 5.33 CeV (b). In (b) the open dots and squares are the 
QCD sum rule predictions of Refs. [O and ITBi, respectively, whereas the full dots 



are the lattice QCD results from Ref. ||T5| 



Fig. 3. The same as in Fig. 2, but for the D Kiu^ (a) and D — > iriui (b) transitions. 
The dot-dashed lines are the monopole approximation with Mpde = 2.11 CeV (a) and 
2.01 CeV (b). The full triangles are the predictions of the lattice QCD simulations of 
Ref. m. 
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Table 1 



weak 


wave 


/+(0) 


/+ (O'max) 


r 


transition 


function 








B ^ D 


GI 


0.765 


1.294 


11.13 




NR 


0.692 


1.322 


9.835 




HO 


0.684 


1.365 


9.780 


B ^ n 


GI 


0.464 


1.803 


15.21 




NR 


0.361 


1.524 


11.12 




HO 


0.293 


1.658 


9.624 


D 


GI 


0.835 


1.379 


0.114 




NR 


0.787 


1.598 


0.111 




HO 


0.780 


1.560 


0.112 


D 


GI 


0.762 


1.281 


0.167 




NR 


0.694 


1.216 


0.156 




HO 


0.681 


1.289 


0.160 



Table 2 



GI 
NR 
HO 



0.464 
0.361 
0.293 



0.762 
0.694 
0.684 



SR |T7| 
SR PI 
SR PI 
WSB H 
ISGW li 
LAT 1T| 
LAT 

Exp PI 



0.29 ±0.01 
0.26 ±0.02 
0.4 ±0.1 

0.33 

0.09 
0.35 ±0.08 
0.43 ±0.02 



0.66 ±0.03 

0.5 ±0.1 
0.75 ±0.05 

0.69 

0.51 
0.80 ±0.08 

0.75l|j:i ± 0.06 
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Table 3 





r(5° n+e-u) 




GI 




1.521 


1.67 


NR 


1.112 


1.56 


HO 




0.962 


1.60 


SR [17] 


0.81 


1.56 


SR [ 


1^1 


0.51 ±0.11 


0.80 ±0.17 


SR 


IS] 


1.45 ±0.59 


1.66l°:i 


WSB 


E 

i 


0.74 


1.41 


ISGW 


0.21 


0.77 


LAT 


m 


0.8 ±0.4 


1 8+0-6 
-•-•0-0.4 


LAT 


|16| 


1.02 ±0.36 




Exp 


111 




-•-•-•-0.6 
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I.L. Grach, I.M. Narodetskii and S. Simula: fig. 1 
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I.L. Grach, I.M. Narodetskii and S. Simula: fig. 
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I.L. Grach, LM. Narodetskii and S. Simula: fi^ 
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I.L. Grach, I.M. Narodetskii and S. Simula: fig. 3b 
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